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ABSTRACT 

The recent transient event Swift J1644+57 has been interpreted as emission from a col- 
limated relativistic jet, powered by the sudden onset of accretion onto a supermassive black 
hole following the tidal disruption of a star Here we model the radio-microwave emission 
as synchrotron radiation produced by the shock interaction between the jet and the gaseous 
circumnuclear medium (CNM). At early times after the onset of the jet (f < 5 - 10 days) a 
reverse shock propagates through and decelerates the ejecta, while at later times the outflow 
approaches the Blandford-McKee self-similar evolution (possibly modified by additional late 
energy injection). The achromatic break in the radio light curve of Swift J1644+57 is nat- 
urally explained as the transition between these phases. We show that the temporal indices 
of the pre- and post-break light curve are consistent with those predicted if the CNM has 
a wind-type radial density profile n oc r"^. The observed synchrotron frequencies and self- 
absorbed flux constrain the fraction of the post-shock thermal energy in relativistic electrons 
e,. ~ 0.03 - 0.1; the CNM density at 10^^ cm nig » 1 - 10 cm"^; and the initial Lorentz factor 
Fj * 10 - 20 and opening angle 0j ~ (0.3 - IWj^ ~ 0.01 - 0.1 of the jet. Radio model- 
ing thus provides robust independent evidence for a narrowly collimated outflow. Extending 
our model to the future evolution of Swift J1644-i-57, we predict that the radio flux at low 
frequencies (v < few GHz) will begin to brighten more rapidly once the characteristic fre- 
quency Vm crosses below the radio band after it decreases below the self-absorption frequency 
on a timescale of months (indeed, such a transition may already have begun). Our results 
demonstrate that relativistic outflows from tidal disruption events provide a unique probe of 
the conditions in distant, previously inactive galactic nuclei, complementing studies of normal 
AGN. 

Key words: black hole physics - galaxies: nuclei. 



1 INTRODUCTION 

Supermassive black holes (SMBHs) are most easily studied when 
they accrete at high rates for extended periods of time, powering 
active galactic nuclei (AGN). Although AGN provide important 
probes of the central regions of distant galaxies, this view is nec- 
essarily biased because periods of bright AGN activity are special 
epochs in the history of a galaxy, during which the substantial ki- 
netic and photon luminosities from the SMBH substantially alter 
the distribution of gas and dust from their more typical "quiescent" 
state. Because the majority of SMBHs are underfed, the conditions 
in most galactic nuclei are thus challenging to study most of the 
time. 
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The high energy transient Swift J164449.3+573451 (here- 
after Swift J 1644-1-57) was detected by the Burst Alert Telescope 
(BAT) on March 28, 2011 as a low SNR image trigger. Sub- 
sequent imaging at radio, optical and. X-ray wavelengths local- 
ized the event to within < 15 pc of the center of a compact 
galaxy at redshift z - 0.353 teerger et alj|201 ll ; iFruchter et al.l 
I2OI ll : llTevan et al]|201 ih . The inferred R-band luminosity places 
an upper limit < 10' on the mass of the SMBH in the galac- 
tic ce nter, if the galaxy follow s the Magorrian et al. (1998) re- 
lation ^Levan et al."2 01ll; iBloom et alji2011,: .Burrows et al.n20lll : 



cf Miller & Gultekinl l201lh . lOuataert & KasenI feOllh point out 
that the SNR of the 5w//j/XRT data is not sufficient to constrain 
significant variability on < 10 s timescales, in which case the mass 
of the central compact object could in principle be much lower, 
possibly consistent with a stellar core collapse event, as are asso- 



©2011 RAS 



2 B. D. Metzger, D. Giannios, P. Mimica 



ciated with normal long-duration gamma-ray bursts (GRBs). How- 
ever, the coincidence of Swift J1644+57 with the galactic nucleus, 
combined with the lack of any previously kno wn GRBs with sim- 
ilar luminosity or duration jLevan et al.ll201 ih . suggest that Swift 
J1644-I-57 most likely originated fr om the rapid onset of accretion 
onto a SMBH jBloom et alj|201 ih . A high energy transient with 
similar properties. Swift J20 58.4-H05 16 (hereaf ter Swift J2058+05), 
was recently discovered by dCenko et al.]2oT ij), who argue that this 
event also resulted from rapid accretion onto a SMBH in a galaxy 

atredshiftz^ 1.18. 

Lack of evidence of prior radio or X-ray activity jLevan et al.l 

l201lh suggest that the rapid increase in the accretion rate responsi- 
ble for Swift J1644+57 and Swift J2058-I-05 cannot result from gas 
entering the sphere of influence of the SMBH (for instance from the 
destruction of a passing molecular cloud), since this would require 
a timescale ~ rsoilcr > 10'* yr to appreciably alter the accretion rate 
near the horizon, where r,o; ~ 1 pc is the radius of the sphere of 
influence and cr ~ 100 km s"' is a typical bulge velocity disper- 
sion. A more plausible source of such a rapid onset of accretion 
is the tidal disruption of a star by a SMBH (tidal disruption event, 
or "TDE": lReeslll98§) . Analytic estimates and numerical calcula- 
tions show that the process of disruption leaves a significant frac- 
tion of the shredded star gravitat i onally bound to t he black hole 
(e.g. iReesI [19881 : lAval et alll2000l ; iLodato et alj|2009 ). The accre- 
tion of this stellar debris has long been predicted to power a thermal 
'flare' at optical, UV, and X-ray wavelengths that lasts for months 
to years after the merg er fe.g. IUlmej[l99^ ;[ Strubbe & Ouataerll 
I2OO9I : iLodato & Rossilbol ih . Until recently, this thermal signa- 
ture was the only means known for detecting TDE candidates 
(e.g. iKomossa & GreineJ 1 19991: iDonlev et alj I2OO2I ; | Gezari et al.l 
I2008ll2009l : lvan Velzen et alj201(]h" 

Just a few months prior to the discovery of Swift J1644-H57 , 
lGiannios' & Metzger (2011) (hereafter GMll) explored the obser- 
vational consequences if a modest fraction of the accretion power 
following a TDE is used to accelerate a coUimated jet to ultra- 
relativistic speeds. Regardless of what ultimately powered the out- 
flows responsible for Swift J1644+57 and Swift J2058+05, if the 
nuclear region was indeed previously inactive then such short-lived 
ejections cannot propagate far from the SMBH before beginning 
to decelerate via their interaction with the dense gas of the cir- 
cumnuclear medium (CNM). This is in stark contrast to long-lived 
AGN jets, which are active for tens of millions of years, dur- 
ing which time the jet evacuates a cavity of kpc scales or larger, 
sometimes depositing its energy far outside of the host galaxy 
(e.g. as in FR II radio sources). GMIl showed that shock inter- 
action between relativistic jet and the CNM powers non-thermal 
radio-infrared synchrotron radiation, which in the case of events 
viewed from a "typical" off-axis angle, peaks on a timescale ~ 
months-years after disruption. GMll did not, however, consider 
the implications of events viewed face-on, as appears to character- 
ize the geo metry of Swi ft J1644+57 (Bloom et al. 201 1) and Swift 
J2058-I-05 fcenkoetal.ll201lh . 

In this paper we extend the basic picture outlined by 
GMll to develop an "afterglow" model for the low frequency 
(radio-microwave) emission from Swift J1644+57. We begin in ^ 
with a basic model for the jet-CNM interaction, using the observed 
prompt high energy emission to estimate the duration and power of 
the jet. We argue that the radio emission from Swift J1644-I-57 is 
synchrotron radiation from behind the forward shock. Using ana- 
lytic expressions verified with numerical simulations, we show that 
the achromatic break observed in the radio light curve is naturally 
explained as the transition between an early phase when the reverse 
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Figure 1. Schematic illustration of the forward-reverse shock structure pro- 
duced by the interaction of a relativistic jet and the sun'ounding CNM dur- 
ing the early phase (1 < ?j) when the reverse shock crosses back through 
the shell of ejecta released during the initial period of jet activity. The for- 
ward shock propagates with a Lorentz factor Ff, into the unshocked CNM 



(of density profile Hci 



), while a reverse shock propagates back into 



the unshocked jet (of density profile «j oc Lj/r^; eq. (3)). A contact discon- 
tinuity separates the shocked CNM from the shocked jet material. Both the 
shocked jet and the shocked CNM move outwards with a similar Lorentz 
factor Fsh , which may be calculated using the relativistic shock jump con- 
ditions assuming constant pressure across the contact (eq. (2))- 



shock crosses the ejecta, to a later phase when the flow settles into 
a self-similar evolution. Using the inferred synchrotron frequen- 
cies and the self-absorbed flux, in 32.1l we constrain the parameters 
of the jet and the CNM. One of our primary conclusions is that 
the radio emission alone provides robust constraints on the Lorentz 
factor and beaming fraction of the jet. In ^we extend our model 
to make future predictions for the radio emission over the com- 
ing months and years. Our one dimensional model is valid until 
the outflow decelerates to mildly relat ivistic speeds, at which poin t 
lateral expansion becomes important IZhang & MacFadvenll2009h . 
We plan to explore this late transition to the Sedov phase, and the 
resulting predictions for the off-axis emission from events similar 
to Swift J1644+57, in greater detail in future work. We summa- 
rize our conclusions in ^ Although our analysis focuses on Swift 
J 1644+57 since more data is available for this event, we conclude 
with a brief discussion of the implications of our model in the con- 
text of Swift J2058+05. 



2 AFTERGLOW MODEL FOR TDE JETS 

We model the interaction between the jet and the CNM guided by 
previous work on GRB afterglows. We adopt a simple model for 
the evolution of the (isotropic equivalent) kinetic power of the jet 



ij.iso(?) = 



hoititj) 



-5/3 



t < tj 
t>t, 



(1) 



motivated by the observed X-ray light curve of Swift J1644+57, 
which shows an initial period of several bright flares of similar 
luminosity, before a se cular declin e consistent with the expecta- 
tion Ljjso oc M oc f ( lReeslll98 8^ from fall-back accretionQ We 
adopt fiducial values for the initial luminosity Lj o = 10''**L48 erg 



' We place hats over the variable t to denote time in the cenfi'al engine (lab) 
frame, thereby preserving the normal variable t for "observer time". 
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Figure 2. Lorentz factor of the shocked circumnuclear medium (CNM) be- 
hind the forward shock (Tsh ; thick solid red line), of the contact discontinu- 
ity {dashed brown line) and of the shocked jet (dotted blue line), extracted 
direct ly from a one-dimensional hydrodynamic simulation iMirnica et al. 
|2009|) . The calculation is performed assuming a jet duration ?j = 3 X 10 
s, initial jet luminosity Lj o = Iff*** ergs s"' (see eq. [T]) and CNM radial 
density profile /icnm — 

10(r/10'**cmr' cm-3 {k = 1). Shown for compar- 
ison with black solid lines are the analytic approximations for Tsh during 
the early phase when the reverse shock is crossing through the shell of 
ejecta released during the initial period of constant jet luminosity (f < tj; 
eq. (2)), and at late times during the Blandford-McKee self-similar evolu- 
tion (; 3> 2/j; eq. (6)) for different assumptions about the total energy of 
the blast wave £j,iso . Note that a break occurs in the Lorentz factor of the 
shocked jet once the outflow reaches the radius r = rcross ~ 3x lO'^ cm (ob- 
server time ? aj fj) at which the reverse shock has crossed entirely through 
the initial shell. The break in Fsh from this transition occurs at a somewhat 
larger radius r » 7 X lO" cm (/ a 2tj), once the rarefaction wave launched 
at r = rc,n"i< reaches the forward shock. 
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Figure 3. Radio- infrared spectra of Swift J1644+57, shown at various 
epochs with approxim ately contemporaneous data. The radio data is from 
IZauderer et al] )201 ih. while the NIR upper limits are L-band detections 
from lLevmetalTj20 llh (also Levan, private communication). Overplot are 
model synchrotron spectra, fit to the data at the earliest and latest epochs, 
where we have assumed that the electron Lorentz factor distribution has a 
power-law index p = 2.3. Notice that both the self-absorption frequency Vsa 
and the peak characteristic frequency v^, evolve to lower values with time, 
in a manner qualitatively consistent with the predictions in equations jlQt 
and mi. 



Figure 4. Radio light curves of Swift J 1644-1-57 from lZauderer et aljj20TT[) 

at several frequencies, most of which are below the self-absorption fre- 
quency at most epochs (cf. Fig. [5). In our model the achromatic break in 
the light curve at ; ~ 10 days occurs when the rarefaction wave (produced 
once the reverse shock crosses the initial shell of ejecta) catches up to the 
forward shock. Shown for comparison are the predicted light curves (arbi- 
traiy normalization) if the CSM density has a constant (dotted line; k = 0) 
or a wind-type (solid line; k = 2) radial profile. 

s"' (see below) and jet duratiorH?^ = lO^/j e s, as constrained from 
the total isotropic X/y-ray energ y £^.150 = 3 x 10^^ ergs measured 
over the first few weeks (e.g. i Levan et"al] I2OI ll : [Burrows et al] 
I2OIII) . Although equation ([T) obviously does not capture the full 
complexity of the observed emission, short-term irregularities in 
the j et power average out in the observed afterglow light curve 
(e.g. Ivan Eerten 



veraee out 
et al]|2009l) 



The inferred expansion rate of the radio emitting plasma from 
Swift J1644-t-57 suggests that the emission originates from a rel- 
ativistic outflow jBloom et aLll201ll : IZauderer et al.ll201lh . while 
energetic arguments suggest that the emission was beamed into a 
bipolar jet. Expanding on the latter point, if the mass of the dis- 
rupted star is ~ 0.5Mq (near the peak of the stellar IMF or similar 
to that of a WD: lKrolik & Piran|201ll) and i f < half of its ma ss ac- 
cretes following the disruption ( IReesll 19881 ; l^al et alj|200oh . then 
reconciling the total inferred energy of the jet x Eyj^e^^ » 
6 X 10^'' ergs with the accreted mass Ma^c ~ Mo/4 requires a beam- 
ing fraction /t = E^°'/E'°^^ ~ 0.01(ej/0.01), where is the true 
jet energy, Cj = £""/A/accC^ is the fraction of the accreted rest mass 
used to power the jet, and ex ~ 0.5 is the radiative effici ency of 
the je t. Motivated by statistics of blazar jets (e.g. Pushkarev et al.l 
12009))^ theoretical cons traints based on models of jet acceleration 
(e.g.l Komissarovlbol ih . and our detailed analysis of the radio data 
below ( 92.1.31 1, we assume that the initial Lorentz factor (LF) of 
the jet obeys Tfi^ < 1, where 6^ is the half opening-angle of the 



^ Note that the jet duration we adopt ~ 3 X lO'' - 10* s is smaller than 
the fiducial value tj ~ lO' s adopted by GMll, who assumed that the jet 
luminosity is constant at times when the fall-back rate is above the Edding- 
ton rate for a ~ IO^Mq BH. GMll also assumed that the pericenter radius 
Rp of the disrupted star was just interior to the tidal disruption radius R,. 
The shorter duration of the bright y-/X-ray flaring from Swift J 1644+57 in- 
stead suggests that either (1) the jet luminosity is directly proportional to the 
fall-back rate, regardless of the accretion Eddington factor, in which case tj 
instead represents the much shorter fall-back time; and/or ( 2) the stellar or- 
bit wa s "deeply plunging" with a pericenter radius i?p <K iCannizzo et alj 
l201lh. 
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jet. In this case the beaming fraction is set entirely by the l/Fj 
emission cone (independent of 6j), such that the value of /b re- 
quired above corresponds to Fj ^ (2/i,)"''^ ~ 10(ej/0.01)"''^. In 
our more precise analysis in 92.11 we determine the initial kinetic 
jet luminosity Ljo (eq. |T|) self-consistently by demanding that 
-Ej.iso = (5/2)Lj_orj = 2£y_i,„(rj6ij)-^ where the factor (Tfi-^fll < 1 
corrects the observed isotropic emission to the isotropic kinetic en- 
ergy of the jet. This correction accounts for the fact that, although 
the jetted emission takes place within a cone of opening angle 
^obs ~ l/Fj, the jet energy is concentrated into the smaller phys- 
ical jet opening angle 9j when ^jFj < 1. 

We now describe the stages in the interaction between the jet 
and CNM using analytic expressions, which we will verify below 
with a full hydrodynamic simulation. Initially after the jet forms, a 
forward shock (FS) is driven into the CNM, while a reverse shock 
(RS) simultaneously propagates back through the ejecta (see Fig.[T] 
for an illustration). Jump conditions across the relativistic shocks 
and pressure balance at the contact discontinuit y give an expressio n 
for the LF of the shocked fluid in the lab frame Isari & Piraril995h : 

rj[i + 2F/" 

where «;.nm is the CNM density and 

I' J 



17 - Ak £j.iso 
16;r m„n,„,^i 



1/2 



(6) 



(2) 



(3) 



is the rest-frame density of the unshocked jet. The LF of the FS is 
given by Ffj = V2Fsi,, while the LF of the shocked material relative 
to the unshocked jet is given by 



(4) 



where the tilde denotes a measurement in the frame of the un- 
shocked jet. 

Equations ([2} and ^ show that the reverse shock is relativistic 
(f,ei » 1) if 2Fj(;7,„,„/nj)-'/2 » i. This condition. 
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(cm-O (lO'^cm)^^^' 



(5) 



is satisfied for reasonable values of the CNM density on radial 
scales probed by the jet during the first few weeks. Below we as- 
sume that the RS is indeed at least mildly relativistic and show that 
this results in a self-consistent model for the radio emission from 
Swift J1644+57. 

Assuming that the RS is (mildly) relativistic, after an observer 
time / ~ tj the RS crosses entirely through the shell of ejecta re- 
leased during the initial period of constant jet luminosity (f < t^). 
After this point a rarefaction wave begins to propagate forward 
from the back of the shell (in effect "communicating" the end of the 
initial period of jet activity to the rest of the outflow). After passing 
through the shocked ejecta and crossing the contact discontinuity, 
the rarefaction reaches the FS on a timescale / ~ 2tj (with the ex- 
act time depending weakly on the speed of the RS). This results 
in a sudden change in the properties of the material behind the FS, 
which we will show below leads to a break in the synchrotron light 
curve. Note that if the RS were not at least mildly relativistic, as 
we have assumed, then the observed break in the light curve break 
would in general occur on timescales much longer than fj. 

Finally, at la te times the outflow approaches the 
iBlandford & McKed ( Il976h (hereafter BM76) self-similar 
evolution, for which the LF of the shocked fluid evolves as 



Note that this expression depends only on properties of the CNM 
and the total ejected energy £j_iso. In the case of TDE, a factor of 
1.5 times more energy is injected at times > fj by the oc r^'' tail 
(eq. U]) than in the initial emission episode {t < tj). Depending on 
when this late material in incorporated into the flow (which depends 
on its velocity with respect to that of the FS), Ej jso will vary from 1 
to 2.5 X fjLj.o. 

Figure |2] provides a comparison between our analytic expres- 
sions for Fsh(r) in equations |2] and [6] with the LF of the material 
behind the FS as determined directly from a one-dimensio nal hy- 
drody namic simulation of the jet-CNM interaction ( Mimica et al.l 
\200% . In this calculation we have assumed a CNM density profile 
n = 10(r/10'^ cm)"' cm"^ (i.e. k = 1) and use the jet power evo- 
lution Lj,iso(0 given in equation l[T) for parameters Lj o = lO"** ergs 
s-',fj =3x10^ s, andFj = 10. 

Figure|2]shows that at early times (f < fj) the LF of the shocked 
CNM (solid red line) from the full simulation agrees well with 
the evolution Ffs oc T^^ oc r"''* predicted by equation Q {solid 
black line) in both slope and normalization. A clear break is ob- 
served in the LF of the shocked jet {dotted blue line) at the radius 
r = ''cros.s ~ 3 X 10" cm at which the RS reaches the back of the 
initial shell of ejecta; as discussed above, this occurs at an observer 
time f » 3 X 10^ s which is similar to the jet duration ?j. A rarefac- 
tion wave is also launched at this time, which reaches the FS at a 
somewhat larger radius ~ 7 x 10'^ cm (observer time ^ 2tj), this 
time producing a break in the LF of the shocked CNM just behind 
the FS, which in turn produces a break in the FS light curve. 

Finally, at times > 2fj, Fsh decreases somewhat more shallow 
with radius than the evolution Ffs oc r"' predicted by the BM76 
self-similar (eq. ||6]), due to the additional energy injection from 
the oc r^^^ tail (which in this example is assumed to have the same 
LF Fj = 10 as the initial jet). Note that the true late-time evolution 
of Fsh(r) lies between the BM76 evolution for a blast wave energy 
2.5Lj_o?j (with oc r^'' tail) and £j = Lj o?j (without tail). However, 
since below we are interested in comparing to data on timescales 
> 2?j, the no-tail BM76 solution (fjiso = i'j.o'j) is used since it 
provides a more accurate estimate of the value of Fsh on timescales 
relatively soon after the break. 



2.1 Modeling the emission from Swift J1644H-57 

Figure |3] shows radio spectra of Swift J1644-F57 at various epochs 
with a pproximately contemporaneous data from IZauderer et al.l 
( I20T 1}), with time measure d starting from the first dete ction of X- 
ray activity on March 26 l iKrimm & Barthelm-^j|201 ih . Note that 
at all epochs the spectrum is self-absorbed Fy oc at the lowest 
frequencies y < 10 - 100 GHz. This indicates that the source is 
relatively compact, consistent with its young age as inferred by the 
lack of evidence for previous activity from the location of Swift 
J1644+57. Shown for comparison in Figure[3]are two model syn- 
chrotron spectra, fit to the earliest and latest epochs of radio data. 
Note that both the characteristic (or peak) frequency Vn, and the 
self-absorption frequency show evidence of decreasing with 
time. Also shown are L-band detections in the near infrared (NIR) 
(with large error bars [not shown] due to possible dust extinction), 
which we treat as upper limits since it is unclear whether the NIR 
emission originates from the same source of synchrotron emission 
as what dominates at radio wavelengths. 

Figure |4] shows the radio light curves of Swift 
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J1644+57 jZauderer et alj l201lh . focusing on several low 
frequency bands which are the least affected by scintillation. The 
data are well fit by a power law Fy oc f with cf « 2 at times t < 10 
days, followed by break to a shallower rise a » 0.5 at later times. 
The apparently achromatic nature of this break suggests that it 
results from a change in the jet dynamics, rather than a spectral 
break passing through a particular bandpass (see also Fig. |3]l. As 
we discuss in the next section, this evolution is consistent with the 
transition expected from the early shell-crossing phase to the late 
BM76 phase around the time t ~ 2tj ~ 10 days, as illustrated in 
Figure|2] Note that for Swift J1644+57 we adopt a value for the jet 
duration ?j « 5 x 10^ s (fj,6 = 0.5). 



2.1.1 Evidence for a wind-type CNM 

Although electrons are (in principle) accelerated at both the PS and 
the RS, in 32.1. 3l we argue that synchrotron emission from the FS 
likely dominates in the case of Swift J1644+57. At frequencies be- 
low the self-absorption break (y < Vj^), the radiated emission is 
given by the usual blackbody emissivity, calculated using a temper- 
ature 3kT ~ y^mgC^ corresponding to that of the minimum LP y,,, 
of shock-accelerated electrons. Here 7ni - (p-2/ p-l)(mp/me)£fT^\^ 
follows from the shock jump conditions, is the fraction of the 
post-shock thermal energy in relativistic electrons, and p a 2 - 3 is 
the power-law index of the electron LF distribution. The resulting 
source-frame specific luminosity is approximately given by 



64;r2 (p - 2) , , , 



(7) 



where the factor (T^^i^dj)^/! < 1 corrects for the missing flux due to 
the finite lateral extent of the jet. 

We now derive how the light curve Ly<^j., scales with time both 
before and after the RS crosses the jet (ignoring the overall flux 
normalization for the moment), in order to constrain the density 
profile of the CNM. If the CNM density scales with radius as a 
power law n^^m °^ '' * (Fig. [Hi, then equations ([2} and ^ show that 
Fsh oc r*"^'''''[oc r<*"^*'-] for the shell-crossing[Blandford-McKee] 
phases, respectively, where we have used equation l[3} and have 
assumed that 2rj -sju^nm/nj » 1 in equation l|2j. Observer time and 
radius are related by f r/2cT^|^, such that r oc ?^'<''"*'[f''*"*"*']. From 
equation Q, we therefore find that Lv<v,, oc r^F^j^ oc ^(t+sj/zjj-tf-i)] 
f<*+2"<*-*'[oc For k = 0(1)[2] we see that the pre and 

post-break temporal indices are given by = 0.5(1)[2] and a2 = 
-l/4(0)[l/2], respectively. As the observed light curves in Figure 
|4] clearly illustrate, the steep initial rise, and continued post-break 
rise, are both broadly consistent with the predictions for a wind- 
type CNM (k ^ 2), but are clearly inconsistent with a constant 
CNM density {k = 0)0 



Note that we have again neglected the possibility of additional energy 
injection from the oc r^/^ tail of the jet luminosity, which may cause F,!, 
to decrease slower with radius (and hence the post-break light curve to rise 
faster with time) than is predicted by the constant-£j BM76 evolution (see 
Fig.|2). However, given the observational uncertainty in the temporal slope 
a2 of the post-break light curve, and the theoretical uncertainty regarding 
when this energy will be incorporated into the forward shock (which de- 
pends on the LF of the jet at late times relative to that of the forward shock), 
this additional complication does not alter our basic conclusions. 



2.1.2 Wind CNM afterglow model 

Assuming that k = 2 with jicnm = '!i8(''/10'* cm)"- cm"^, equations 
(|2) and (|6) can now be written 



1/4-1/4 



l/2,l/2_^-l/2_-l/2 

"is 



10.9L;-r;.-«|--r 



: 7.2L 



1/4 -1/4 Z^) 
48 "18 j 



t < tj 



t > 2tj 



(8) 



where ntg is in cm and in the last expression we again relate the 
emission radius r = rig 10'** cm and observer time by the expression 
r^ 2rl^ct/(i +z), viz. 



^8 



0.1s(l)t,,L\i^n-l'^ 
2.3 (t) 



t<ti 



, 1/2 -1/2 , ^ 
,6^48 «ig' t > 2tj 



(9) 



Again note that we have assumed £j = Ljtj (i.e. neglecting the 
tail) in the second expression in equation l[8j since we are in- 



oc r5/3 

terested in the jet properties on timescales of < few xfj (see Fig.|2]l. 

Based on the evolution for T^hit) in equation we now give 
expressions for the time evolution of the characteristic synchrotron 
frequencies and the self-absorbed flux. First, the peak (or charac- 
teristic) synchrotron frequency is given by 



2TmicC{\ + z) 

6xlO'OHz<_,4/^,Liffr^(i)"', 

2-^48 fj.6 j 



1.7 X 10" Hze; ,e, 



1/2 



-3/2 



t < tj 



t > 2tj 



(10) 



where B = s/Sne^i^ and e,h - AT^^^nipncnmC^ are the magnetic field 
and energy density of the shocked fluid, respectively; = 0. le^ -i; 
ee = 0.01eB,-2; and we have assumed an electron LF power-law 
index p = 2.30 

Likewis e, the tim e evolution of self-absorption frequency is 
given by fe.g. lOranot & Sari 2002) 



3x 10''Hze-i,4''2.^48''^.^ 



-1 1/5 r-2/5,- 
•^f.-l "^8,-2^48 



3/5 



t<ti 



t > 2tj 



(11) 



Equation Jilt assumes that Vsa < Vm, which Figure[3] shows is valid 
for Swift J1644-I-57 at times ? < 23 days (In particular the spectral 
slope above v^^ a; 10 - 100 GHz is consistent with Fy oc v''^, but is 
much different than Fy oc y^'-, as would be expected if Vn, < Vsa). 

Finally, the flux of the self-absorbed emission at the distance 
of Swift JI644+57 (z ^ 0.35; = 5.7 x 10" cm) is given by 



K(v < Vsa)(l +z)- 
4nDl 

26 mJye^o.ii^ 



24OmJye,,o.,L^fr^^n78^'^v?„0j_,(l) 



t < ti 



(12) 



t > 2tj 



where 9j = 0.19 v = lO'^vio Hz, and we have corrected the nor- 
malization of the flux in second expression to account for the de- 
ta iled radial structure of the outflow, as described in the Appendix 
of lGranot&SMil ( l2002h . 

A casual comparison of the above expressions to the evolution 



Although we assume p = 2.3 throughout this paper, the value of p always 
appeal's in the combination jj^rj] fc. such that one can adopt a different value 
of p thi'ough a simple renormaUzation of the value of e^. 
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n,8 (cm"^) n,8 (cm"^) 



Figure 5. Constraints on the fraction of tlie shock energy placed into rel- 
ativistic electrons (eg) and the CNM density at r = lO'** cm (nig) from 
fitting the synchrotron spectrum from the final epoch of radio observa- 
tions on April 17 {t x 4/j) to our FS model. Constraints shown include the 
peak frequency (blue region; eq. 1101 ). the self-absorption frequency Vja 
(brown region; eq. 1111 ). and the self-absorbed flux at 4.9 GHz Fv=4.9GHz 
{black line; eq. 1121 ). including estimates of the uncertainties in the mea- 
sured quantities. In this example we take characteristic values for the frac- 
tion of the shock energy in magnetic fields eg, = 0.01 and the parameter 
= Fj^j = 1, which results in a unique solution {orange star) for the val- 
ues of ff, X 0.15, «i8 X 1.5, and the initial LF of the jet Fj = 17 (see also 
eqs. I14I16I ). The green fine shows the upper limit on the density, such that 
the FS dominates over the RS emission as obtained using equation )19t . 



of the model synchrotron spectra in Figure [3] shows good overall 
agreement. For example, is observed to decrease from > 3x 10" 
Hz at f 5 days ~ fj to ~ 3x 10'° Hz at f 23 days x 4fj , consistent 
with the prediction Vn, oc r^'^ for t > 2/j from equation dlOt . Like- 
wise, Vsa decreases from ~ fewxlO'" Hz to ~ 10'" Hz, consistent 
with the slower evolution v^.^ oc r^^^ predicted by equation Jl It . 

We now use the synchrotron frequencies and self-absorbed 
flux from the last epoch of data on April 17 (/ a 23 d a 4?j) to 
constrain the unknown parameters 6^,^18, and Fj. As discussed near 
the beginning of this section, the observed high energy fluence con- 
strains the isotropic jet luminosity to obey 



0.577-2, 77 < 1 



0.5, 



77 > 1 



(13) 



where 77 = Tfiy Throughout this section we assume that 77 < 1, 
although in 32.1.3| we consider (and argue against) the possibility 
that 77 > 1 . 

Adopting fiducial values for the characteristic frequency 
y„,(f = 23 d) = Vm,23 ~ 40 GHz, self-absorption frequency Vsa(? = 
23 d) = y,a,23 ~ 10 GHz and the 4.9 GHz flux F^=4.9ghz(? = 23 
d)= F4.9,23 ~ 2 mjy motivated using the data in Figure[3] equations 
lll0t-fT2t can be solved for the following 



: 0.126: 



ni8 = 1.5e, 



"MlOGHz/ ' 
-'UOGHz/ \ 10 GHz/ ' 
' B-n40GHz/ llOGHz/ 



2mJy 



(14) 
(15) 

(16) 



10 



Equations lll4t-lll6|l show that for rea sonable values of « 
^ - 0.1 (e.g. lPanaitescu & Kumaill200lh and 77 = FjSj ~ 0.1 - 1 



Figure 6. Same as Figure|5] but calculated for different values of 6b = 10 
and;; = Yfi^ = 0.1. 



(e.g. IPushkarev et 'Zl l2009h . one finds reasonable values for eg ~ 
0.05-0.2, 7118 ~ 0.3-10, and Fj ^ 10-20. Note that the value of Fj 
is particularly robust since it depends only weakly on the uncertain 
values of 6b and 77. As we discuss further in ^ this implies that 
the beaming fraction of Swift J 1644-1-57 is well constrained to be 
/b ~ l/2r^ « 1 - 5 X 10-3. 

Figures [5] and |6] show the parameter space of and 7ii8 con- 
strained by the observed values of yni,23, v'sa,23, and F4 C)_23, includ- 
ing estimates of the uncertainties in these quantities. We show two 
cases corresponding to different values of and Fj, the latter cho- 
sen from equation M6\ to produce an acceptable solution for which 
the constraints on y^a, ym, and F,,=49ghz all overlap. 



2.1.3 Constraints on the value ofrj = FjSj 

Although Fj is well-constrained by equation 1 II6I 1, the combination 
77 = FjSj is poorly constrained by the analysis thus far. In this section 
we present several independent arguments for why the value of 77 is 
~ 1. 

We begin by placing a lower limit on 77. One constraint not 
yet considered is that for the RS to be at least mildly relativistic, 
as is necessary to explain the timing of the observed radio break, 
then the second term in parenthesis of Equation l[2j must be » 1 
(see eq. fS]). Requiring that Fsh <K Fj at times t < tj results in the 
following constraint: 



,-,r,>o.oH-(j^)-"'"( 



Vsa.23 \ 

10 GHz/ 



10/9 / 



I -^4.9.23 \ 



4/3 



,(17) 



\ 2 mJy / 

where we have used equations il5i and J16t . 

Although the constraint in equation l ll7t is relatively weak, 
one can motivate a somewhat higher value for 77 by considering the 
physical value of CNM density encountered by the shock at the last 
epoch of observations. Using equation l|9j we find that 



: 0.477^,772 



4fi 



(18) 



which shows that if 77 ~ 0.1 (implying 7!i8 < 6 for 6b > 10"' using 
eq. 1151 ), then the density probed at f a 4fj is 7!c„,n < 0.1 cm"'. If 
this value is accurate, then the CNM density in the host of Swift 
J1644-1-57 is much lower than the density > 10 cm '^ inferred on a 
similar radial scale in the center of the Milky Way jBaganoff^ et al.l 
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l2003l : IOuataertl2004[) . Such a low density would be particularly sur- 
prising because several magnitudes of dust extinction were inf erred 
along the line of sight to Swift J1644+57 l lBloom et al1l201ll) . We 
conclude that the unphysically low CNM density implied by ;y « 1 
instead favors a value for r] >03. 

An upper limit on the value of 77 can be placed by requiring 
that emission from behind the FS dominates that from the RS. The 
light curves of both the FS and the RS at frequencies v < Vsa ini- 
tially rise as Fy oc f for a wind-type CNM, such that the data at 
times f < fj cannot alone be used to determine which dominates 
the emission. After the shock crosses the ejecta (/ > /j), however, 
the RS emission is predicted to decrease more rapidly (or increase 
less rapidly) than is observed in Swift J1644-I-57 . This is because 
no new electrons are accelerated at the RS at times t > tj, yet the 
plasma loses thermal energy to adiabatic expansion. In particular, 
emission fro m the RS at v > Vm is predicted to decrease as Fy oc 
or faster (e.g. lSari & Piranll999t) . contrary to the ob served high fre- 
quenc y post-break light curve of Swift J1644-I-57 jZauderer et all 
l20Tlh . 

On timescales t ~ tj the brightness of the forward shock rel- 
ative to that of the RS is given by the ratio of the relative electron 
temperatures (see eq. (7)), which is in turn proportional to the ratio 
of the relative LFs 



Wind CNM, t^,„(25 d)=10'° Hz, i/^(25 d) = 5x10'° Hz 



Ly.yJFS) 2Tl 



~ 2.5n,. 



(19) 



where we have used equations © and {8]( for L48 = 0.5 (eq. 1131 ). 
Again solving for the unknowns as in equations ( 114116b above, but 
this time assuming 77 = Fjffj > 1, we now find 



: 0.126: 



1 30 GHz/ 

, . -3/8/ Vm.23 V^'V na,23 

,^ = o.o5.-f f^:^)'"V^:^r 



^4.9 



2mJy 



1/2 



(20) 
(21) 
(22) 



Substituting the expressions for njg and 6j into equation l ll9t gives 

Ly.vJFS) 



Ly^yJRS) 



1 3 ,/32 / JW^r''^« ^ JV2^f IF^V" 

B,-n40GHzj llOGHz/ UmJyj ^" 



For the 

'7 = m 
To 

range ~ 
density 
of r] are 
jets (e.g. 
models 



FS to dominate over the RS (as observed) we thus require 
< 1. 

summarize, the available data constrain 77 = Fjffj to lie the 
0.05 - 1, but given a realistic lower limit on the CNM 
77 most likely lies in the range ~ 0.3 - 1. These values 
consistent with those inferred from observations of blazar 
JPushkarev et alj2009l) and those expected from theo retical 
of relativistic jet acceleration (e.g. lKomissarovll201 ih . 



3 FUTURE EVOLUTION 

3.1 Ongoing Blandford-McKee Phase 

In this section we extend our model for the synchrotron emission 
from Swift J1644+57 to predict the ongoing evolution of the ra- 
dio light curves over the coming months to years. If the outflow 
continues to evolve as predicted by the BM76 evolution with con- 
stant energy in a wind-type CNM (eq. ||6]), then the self-absorbed 




100 

Time since March 26 (days) 

Figure 7. Future evolution of the radio light curves of Swift J 1644+57, 
assuming that Vja = lO'" Hz and Vn, = 5 X lO'" Hz at f = 23 days, and that 
the CNM continues to maintain a wind-type radial density profile 77i;nm 



Wind CNM, i's<,(23 d)=10'° Hz, !/^(23 d) = 2x10'° Hz 



"^'^■^10 CHz 



012 ^ 



100 

Time since March 26 (days) 

Figure 8. Same as Figure [T] but assuming that y^, = 2 X lO'" Hz at f = 
23 days. Taking into account uncertainty induced by scintillation at high 
frequencies, note that the flux at 90 GHz does not match precisely onto its 
observed value at ( = 23 days in this case. 

Const CNM, t^,„(25 d)=10'° Hz, i^^(25 d) = 2x10'° Hz 



0.1 



90 GHz 



012^ 



100 

Time since March 26 (days) 

Figure 9. Same as Figure [8] but instead assuming that the CNM density 
profile flattens to a constant value at times t ^23 days. 
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flux and peak frequency will continue to evolve as Fy^y^^ oc f''^ 
(eq. 1121 ) and Vn, oc r''^ (eq. 1101 ). respectively. The self absorption 
frequency will continue to evolve as v^a « f^'^ while < Vm, but 
then will chan ge its temporal slope slightly to v^a oc f-(3p+5)/(2p+io) 
once Vsa > Vm lOranot & Saril2002h . 

Using this information alone, Figures |7] and [8] shows predic- 
tions for the future evolution of the radio flux at several common 
observing frequencies, for two different assumptions about the peak 
synchrotron frequency at f = 23 days (which, as shown in Fig. [3] is 
somewhat uncertain). At high frequencies > 10 GHz, the flux will 
continue to decrease monotonically with time, approaching < 0.1 
mjy at 90 GHz on timescales of a year. At low frequencies the 
evolution is more interesting: at 1.4 and 4.9 GHz the flux should 
initially continue to rise with Fy oc r''^. Then, on a timescale of 
months, the peak synchrotron frequency will decrease below the 
self-absorption frequency, after which point the spectrum is rises 
as Fy oc for Vm < v < Vsa. When an observed bandpass falls be- 
tween this frequency range, the flux is predicted to rise as Fy oc t^^* . 
Since this condition is likely to be satisfied at 1.4 GHz (and possi- 
bly at 4.9 GHz), the temporal slope of the light curve in these wave- 
bands may suddenly increase by a value > 1/2 after a timescale of 
months. 

Although the CNM density inferred from the early radio data 
is consistent with a wind-type radial profile, the CNM profile could 
in principle change with radius, in particular considering that very 
low CNM densities are implied by the wind-like profile at late 
times (see eq. 1181 and surrounding discussion). In order to explore 
this possibility, we also consider how the radio evolution would 
change if the CNM density flattens to a constant value at times 
> 4?j X 23 days. In this case, the evolution can be calculated in 
the same manner as above, but using scalings for the synchrotron 
flux and frequencies which are instead appropriate for a constant 
CNM (k = 0): Fy^y^.^ oc y,„ oc r^/^ v^a oc (when v,., < y„,). 



cx t 



-(3p+l)/(2,;+10) 



(when Vsa > Vm)- Figure |9] shows our results 
in the constant CNM case, which show that initially the flux at all 
frequencies begins to decrease. However, this time when the self- 
absorption frequency crosses below the characteristic frequency, a 
low frequency rebrightening occurs on a timescale of months. 

Thus far we have neglected the effects of the oc r^^^ tail of the 
jet luminosity (eq. |T|) on the late-time radio emission. Depending 
on when this material catches up to the forward shock, energy de- 
position from the tail may cause Fsi, to decrease less rapidly with 
time than would be predicted with no additional energy injection 
(see Fig.|2j. The effect of the tail in this case may be to increase the 
temporal indices of the light curves above those predicted in Fig- 
ures |7l9] Though more speculative, another possible source of late 
emission may occur from the RS, if it continue to propagate back 
through t he ejecta in the t ail after passin g through the initi al ejecta 
shell (see lUhm & Belobor odov 20071 and lOenet et al.l2o"07l for such 
a model applied to normal GRB afterglows). A detailed exploration 
of this idea is, however, beyond the scope of this work. 

3.2 Sedov Phase 

Our calculation thus far assumed an ultra-relativistic outflow, which 
allowed us the simplification of treating the flow as essentially 
one-dimensional. Recent 2D hydrodynamical simulations of de- 
celeration of a collimated blastwave demonstrate that the flow re- 
mai ns approximately conical un til t he LF of the FS de creases to 
< 2 dZhang & MacFadverj|2009l : cf. lWygoda eTalllloT 11) at a time 
t ~ fNR- This conclusion is supported b y recent analytic mod- 
els ( lLvutikovll201 ll ; ICranot & Piranll201 ll) . After a transient phase 



of lateral expansion on timescales t > ?nr, the blast will ulti- 
mately relax into a non-relativistic, spherical expansion (Sedov- 
Taylor phase) centered at the location of the deceleration of the 
blast (GMU). 

Using equation ^ (now increased by a factor 2.5''^ ~ 1.6 to 
now account for the additional energy injected by the oc r^^^ tail) 
we estimate that the non-relativistic transition (taken to be F^h g 2) 
will occur at a time 



fNR(A: = 2) - 8.8yr;7-2„7i « 5.9 yre^^'^'?"''"' 



(24) 



where in the second expression we have substitute equation il5l for 
nig. Depending on the value of r], the Sedov phase is thus predicted 
to begin on a timescales of years to a decades. 

The non-relativistic transition may actually occur somewhat 
earlier than predicted by equation l l24t if the density of the CNM 
does not decrease oc r^^ all the way out to the Sedov radius, as 
would be expected if the ambient ISM places a density "floor" 
«cnm <; Wism ~ 0. 1 - 1 Cm"'. Substituting ?nr into equation il8t . 
we infer that the physical density at the Sedov radius is ncnm('NR) ~ 
3x 10^'njg cm"' for r] ~ 1. Thus, depending on the values of «i8 and 
Wisra, the non-relativistic transition may instead at the usual Sedov 
time for a constant ambient density: 

t^^(k = 0)^3.6yru-"'(^Y'\ (25) 
\ cm"^ / 

somewhat earlier than predicted by equation l l24t 

The angular size of the radio emission at the non-relativistic 
transition can be estimated to be 



2e,r(t = ?nr) 
= X 0.166j-i«ig mas. 



(26) 



where Z)a ~ 1-0 Gpc is the angular diameter distance to Swift 
J 1644+57 . Here we have used equation l9ll to estimate the jet ra- 
dius at the non-relativistic transition r(t = /nr), assuming (opti- 
mistically) that the latter occurs at the time given by equation l l24t . 
Unfortunately, even for optimistic assumptions, the angular scale 
implied by i26l is unlikely to be resolvable with VLBI, especially 
given the low predicted flux at high frequencies at late times (Fig.jTJ 



4 CONCLUSIONS 

Building on the work of GMll, in this paper we have developed 
a model for the radio emission from Swift J 1644-1-57 produced by 
the shock interaction between the jet and the gaseous CNM. We 
have shown that the achromatic break in the observed light curve 
(Fig. |4) is well explained as the transition from the early reverse 
shock crossing phase to the late BM76 evolution (Fig. |2]l if the 
CNM has a wind-like radial profile (Fig. |4ll. One implication of 
a wind-like CNM is that the post-shock LF is constant during the 
shell-crossing phase (/ < fj) and only decreases slowly with time 
Fsi, oc r'/"* thereafter (eq. (8)). The approximately fixed value of 
Fsi, ~ few that we find over the first few weeks is one explana- 
tion for the approximately constant expansion velocity inferred by 
Zauderer et al. (2011). We emphasize, however, that the initial LF 
of the unshocked jet Fj is much higher than its post-shock value. In- 
deed, a value of Fj » I is required for beaming of the initial X-ray 
emission, as is necessary to reconcile the isotropic energy of Swift 
J1644+57 with the accretion of a stellar mass object. 

Achromatic breaks in afterglow light curves can in principle 
also occur when the LF of the shocked jet Fsh decreases below ~ 
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1 /9j (a so-called "jet break"), as is commonly discussed in the con- 
text of GRBs (e.g. Frail et al. 2001). A standard jet break cannot, 
however, explain the radio break observed in Swift J1644-t-57 for 
three reasons. First, as discussed in 32.1.31 the absence of domi- 
nant emission from the RS constrains the value of = Fj^j to be 
< 1 (eq. 1231 ). in which case Tg^dj is < 1 even at early times. Sec- 
ond, for a wind-type CNM (as inferred from the pre-break light 
curve), no change in the temporal slope is predicted after a putative 
jet break because Fji, is constant with time during the RS cross- 
ing phase. A jet break that occurs after the RS crossing is ruled out 
because it would then be the second break and only one break is ob- 
served. Finally, if the observed radio break were in fact a jet break, 
then its occurrence at the time fa; 10 d a; 2/j would be entirely co- 
incidental. In our model the observed time of the break is naturally 
explained as when the rarefaction wave reaches the forward shock 
(Fig.O. 

Using the inferred synchrotron frequencies and the self- 
absorbed flux, our analysis results in reasonable values for the 
fraction of the post-shock thermal energy in relativistic electrons 
e, ~ 0.05 - 0.2; the CNM density at 10"* cm 1 - 10 cm^^; 

and the initial LF of the jet Fj ~ 10 - 20 (eq. (Hl-flll). Although 
the jet opening angle is less well-constrained, it most likely lies in 
the range ej ~ (0.3-l)xrri _ o.Ol-O.l (^2X3}, with much lower 
values disfavored since this requires an unphysically low value of 
the CNM density (eq. 1181 ). Although our model makes several 
simplifying assumptions (a jet with a sharp edge; broken power- 
law spectra; constant microphysical parameters), it is nevertheless 
gratifying that the properties of the relativistic outflow that we infer 
are remarkably consistent with those thought to characterize nor- 
mal blazar and AGN jets. 

Perhaps our most important conclusion is that the radio emis- 
sion from Swift J 1644-1-57 provides evidence for a narrowly colli- 
mated jet, which is independent of (and stronger than) other argu- 
ments based on rates or total energetics. In particular, the beaming 
fraction of the jet fi, = 6^12 where 9b = niax[Fj"', Sj] is constrained 
to lie in the relatively narrow range ~ 1 - 5 x 10"^ (see eqs. 1161 and 
1221 ). This conclusion is robust and independent of the assumed 
value for r] = Tfiy 

The beaming fraction that we infer for the jetted emission 
from Swift J1644+57 has several implications. First, it implies 
that the true beaming-corrected peak luminosity of the prompt X- 
ray/y-ray emission may be as low as ~ lO"*"" erg s"', similar to the 
Eddington luminosity of a 10' Mq SMBH. Super-Eddington accre- 
tion may thus not be required to explain the observed X-ray emis- 
sion, whether it originates internal to the jet or from the disk itself 
jSoc rates 2011). On the other hand, the beaming-corrected energy 
of the event was ~ 10" ergs, which in the case of jetted emission re- 
quires a jet efficiency of 6, a; 10"^ if the event was indeed powered 
by the accretion of a solar-mass star. Note also that the beaming 
fraction we infer appears to be larger than that required to reconcile 
the rate of Swift J1644-l-57-like events with the TDE ra te inferre d 
by independent me ans (e.g. Wa rig & Merrittl2004 l van Velzen et al.l 
1201 ll ; lBoweil201 ll;[kesden.,201 1 ) by a factor of ~ 1 - 30 (see dis- 
cussions in lBurrows et alj|201ll Cenko et"aill201ll) . This suggests 
that only a small fraction of TDE events are accompanied by [de- 
tectable] relativistic outflows. Finally, the narrow beaming angle 
that we infer may in principle also constrain the spin of the SMBH 
via the [lack of] spin-induced precession inferred by the approxi- 
mately steady jet luminosity ove r the first several m onths of obser- 
vations jStone &Loe"bll201ll ; cf. lLei & Zhan3l201ll) . 

In ^ we make predictions for the future evolution of the 
radio-microwave light curve of Swift J1644H-57 over the com- 



ing months-years (Fig. l7l9l >. At high frequencies the radio emis- 
sion will continue to decrease with time, while at low frequencies 
the predicted light curves are initially relatively flat, depending on 
whether the CNM density continues to be well-described by a oc r"^ 
radial profile (or whether the profile flattens) and when/whether ad- 
ditional energy reaches the FS from the oc r^'^ tail. Interestingly, at 
the lowest frequencies < few GHz we predict that the flux may be- 
gin to increase more rapidly (or even rebrighten) on a timescale of 
months, once the characteristic frequency decreases below the self- 
absorption frequency (indeed, such a transition may already have 
begun). These predictions can be used to distinguish our afterglow 
model for the radio emission powered by external CNM interaction, 
from other models whi ch invoke synchrotron emission entirely in- 
ternal to the jet itself jvan Velzen et alj|201 Tl; iMiller & GultekinI 
[2011). 

An important theme of this work is that relativistic outflows 
from tidal disruption events provide a unique probe of the condi- 
tions in distant, inactive galactic nuclei. Such observations are rare, 
except in the case Sgr A* and nearby low-luminosity AGN, and are 
unprecedented in galaxies at high redshift. If the wind-type CNM 
that we infer from Swift J1644-I-57 truly indicates an outflow from 
the nucleus of the host galaxy, then the mass loss rate [prior to the 
TDE] may be estimated as 



M = Anrv^n^^^mp - 3 X lO^^Mo yr-'^jg (-7--^^^=^ r) , (27) 

\300kms"'/ 

where v„ is the wind velocity, scaled to a value characteristic of the 
expected escape speed from the SMBH at radii ~ 10'^ cm probed 
by the radio afterglow. Depending on the value of «i8 ~ 1 - 10, 
this mass loss rate is similar to that expected from ~ 1 - 10 Wolf- 
Rayet (WR) stars. If scaled to the < 10' total number of WR stars 
in the Milky Way, this would suggest that the steady-state star for- 
mation rate is ~ 10"' - lO^-M© yr"' within the inner pc (given 
the Milky Way SFR rate ~ Mq yr"'). Alternatively, our inference 
that the mass loss rate is si milar to that of a massive star is also 
consistent with the model of 'O uataert & Kasenl ( l20Ill) . who argue 
that Swift J1644-)-57 may simply be an unusual long-duration GRB 
associated with a stellar core collapse event. 

Finally, given the uncertainties, the radio data may also be 
consistent with an ncnm °^ ''^'^ density profile, as predicted for 
Bondi accretion. In this case the accretion rate that one infers is 
similar to the outflow rate given in equation Ml\ . Interpreted this 
way, M is several orders of magnitude lower than the accretion rate 
required to explain the outburst from Swift I1644-I-57. This would 
provide additional evidence that Swift I1644-I-57 resulted from a 
very large increase in the SMBH accretion rate over its long-term 
average, as likely can only be explained by the tidal disruption of a 
star. 

Although we have focused on modeling the emission from 
Swift J1644+57, our results can also be app lied to the other re - 
cent jetted TDE candidate Swift I2058-I-05 jCenko etal]|201lb . 
The X-ray and radio brightness of Swift J2058+05 is similar to 
that of Swift J 1644+57 at a similar epoch, except that no spectral 
steepening is observed at low frequencies, suggesting that the self- 
absorption frequency obeys < 3 x 10' Hz. Equation l ll5t shows 
that one way this observation could be explained in our model is if 
the CNM density is lower by a factor of a few than in the case of 
Swift J 1644+57, possibly consistent with the lower optical extinc- 
tion inferred towards Swift 12058+05. 
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